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Abstract 

Broadening of the transverse momentum of a parton propagating through a medium 
is treated using the color dipole formalism, which has the advantage of being a well 
developed phenomenology in deep-inelastic scattering and soft processes. Within this 
approach, nuclear broadening should be treated as color filtering, i.e. absorption of 
large-size dipoles leading to diminishing (enlarged) transverse separation (momentum) . 
We also present a more intuitive derivation based on the classic scattering theory of 
Moliere. This derivation helps to understand the origin of the dipole cross section, 
part of which comes from attenuation of the quark, while another part is due to mul- 
tiple interactions of the quark. It also demonstrates that the lowest-order rescattering 
term provides an A-dependence very different from the generally accepted A 1 / 3 behav- 
ior. The effect of broadening increases with energy, and we evaluate it using different 
phenomenological models for the unintegrated gluon density. Although the process is 
dominated by soft interactions, the phenomenology we use is tested using hadronic 
cross section data. 
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1 Introduction 



A high-energy parton propagating through a medium experiences multiple interactions that 
increase its transverse momentum. This broadening of the transverse momentum can be 
measured in different reactions, a few examples of which follow. 

In the Drell-Yan process []TJ the lepton pair produced carries undisturbed information 
about the transverse momentum of the projectile quark which undergoes initial state inter- 
actions. The important condition is the shortness of the coherence time of the Drell-Yan 
process @, which allows one to factorize the initial state interaction from the cross section 
for the Drell-Yan reaction. In the opposite regime of very long coherence time, the effect 
of transverse momentum broadening also exists, although the lepton pair is produced as a 
fluctuation long in advance the nucleus (in the nucleus rest frame, see @) and does not 
"communicate" with the quark any more. Nevertheless, the nucleus supplies the fluctuation 
with a larger mean value of momentum transfer than a nucleon target, therefore it is able 
to liberate harder fluctuations, i.e. those which have larger intrinsic transverse momenta 
fcr- As a result, the value of (k^) of the lepton pair increases. The nuclear modification of 
the transverse momentum distribution of lepton pairs in the limit of long coherence time is 
calculated in || £|. 

In a similar way, the production of heavy quarkonia off nuclei QXJ] can measure the 
broadening of a projectile gluon. Final state elastic rescattering of the produced quarkonium 
can be neglected since the cross section is very small. 

One can also study the broadening of a quark originating from DIS on a nuclear 
target, which also includes the Fermi motion of the participating nucleon. Although the 
available data from the EMC experiment || do not show any significant effect, this is related 
to a suppression by a factor z 2 , where z is the fraction of the quark momentum carried by 
the produced hadron. In this experiment (z) ~ 0.25. Nevertheless, new precise data are 
expected soon from the HERMES experiment. 
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In the production of two high-pr back-to-back hadrons (or jets) off nuclei, one of the 
hadrons defines the scattering plane, while the acoplanarity of the other one serves as a 
measure for the nuclear broadening of the transverse momenta by initial and final state 
interactions. Available data || [?| demonstrate an unusually strong effect (see the interpre- 
tation in 0). 

The propagation of a high-energy parton through a medium can be treated intuitively 
as a random walk in transverse momentum space leading to a linear increase of (k^) with 
the thickness of the matter covered ||, 10]. Such an interpretation faces, however, certain 
difficulties. The cross section for the interaction of a colored quark diverges at small mo- 
mentum transfer, and one has to introduce an infrared cut off. On the other, hand the mean 
value of the momentum transfer squared diverges at the ultraviolet limit. 

It is demonstrated in Section |2| that the broadening of the transverse momentum of a 
parton propagating through a nucleus can be treated within the light-cone dipole approach 



introduced in ]TT[ as a color filtering effect ][E| for a color qq dipole traveling through nu- 
clear matter. Color filtering is a size-dependent attenuation of qq color dipoles propagating 
through nuclear matter. 

Since the formal transition from a single quark to a qq dipole propagating in a medium 
might be difficult to understand, a more intuitive derivation based on the classic multiple 
scattering theory of Moliere is presented in Section |3|. It goes along with the treatment 
of multiple interactions in the abelian case by Levin and Ryskin ||13[] . Although it does 
not have explicit QCD input, the final result has the form of a color-screened dipole cross 
section, and one can see where different parts of this cross section come from. One can also 
trace the origin of the /l 1//3 -dependence of {k 2 -) broadening to multiple interactions and see 
that the lowest-order contribution actually has a quite different behavior. 

The energy-dependent dipole cross section is evaluated in Section [| via the unintegrated 
gluon density, which is estimated within models adjusted to describe data for the proton 
structure function over a wide range of x and Q 2 , as well as for the hadronic total cross 
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sections. The predicted broadening is expected to rise steeply with energy. 

Since fcy broadening of a high energy parton results from multiple gluonic exchanges with 
nucleons, nuclear shadowing of gluons diminishes the effect of broadening. This can be also 
interpreted as the Landau- Pomeranchuk effect, or coherence-suppressing gluon radiation, 
which gives an important contribution to the broadening. In Section [| we evaluate the effect 
of shadowing relying on the light-cone Green function formalism, which includes the strong 
nonperturbative interaction of gluons that dramatically diminishes the effect of shadowing. 
The reduction of broadening due to gluon shadowing turns out to be rather small, less than 
10%. 

Nuclear broadening of kx can be affected also by the finiteness of the experimental 
aperture. In Section ||we estimate the corresponding correction factor, which is rather close 
to unity if the aperture covers a range of exceeding a few GeV. 

The main observations are summarized and discussed in the last section. 



2 Transverse momentum broadening as color filtering 

Nuclear broadening of the transverse momentum of a quark propagating through a medium 



was treated in [TJ| in terms of eikonalized multiple interactions of a colorless qq pair via 



gluon exchanges. It was found that the mean transverse momentum squared grows as 

5(k 2 T ) = 2C p A L . (1) 

Here L is the length of the path covered by the quark in nuclear matter up to the point that 
the lepton pair is produced, and pa is the nuclear density. The factor C originates from the 
expression for the total cross section for the interaction between a nucleon and a colorless 
qq dipole having transverse separation tt and cm. energy squared s fLT|l , 



0qq(rT, s) = C(r t s)r 2 T . (2) 
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The r\ behavior at small ry — > is dictated by gauge invariance and the non-abelian nature 
of QCD. The energy-independent part of C was estimated in the Born approximation to 



be C ~ 3 [IT , [15]. However higher order perturbative corrections lead to a rising energy 
dependence of the C(rx, s). At small tt the factor C(tt, s) is related to the gluon density 
in the proton p]j , 

C(r T ,s) = ^G(x,Q 2 ) , (3) 

where G(x, Q 2 ) = x g(x, Q 2 ) is the gluon distribution function which depends on Q 2 ~ l/ r T 
and x = Q 2 /s. The the broadening corresponding to (|I]) takes the form, 

S(k 2 T ) = ^G(x,Q 2 )p A L . (4) 



This expression is also derived in fll~8 



Although the result ([lj) is correct it was poorly motivated in |TJ]]. Better derivations 



one can find in JT8|, [L9|. They, however, assume Gaussian shape for the kx distribution 
and employ the approximation of constant C(tt) which cannot be justified for soft multi- 
ple interactions. Here we derive a general expression fll~6|) for nuclear modification of the 
transverse momentum distribution which is free of these approximations. 

The effect of the mean-square transverse momentum of the quark after propagation 
through a medium depends on the reaction, nevertheless, the broadening S(k^) is universal. 
This is demonstrated below for the example of a hadron-nucleus interaction where we are 
interested in the final transverse momentum distribution of one of the projectile quarks. 
This can be expressed in terms of the density matrix of the final quark, b'), where b is 
an impact parameter, 

dNq _ f 727 727 /. 



J d 2 b d 2 b ! exp [i k T (6 - b')] Q}{b, b') . (5) 



d 2 k T 

The distribution dN q /d 2 kT is normalized to one. The final density matrix is related to the 
initial one as 

Wb,b') = TrS\b' + B)nl(b,b')S(b + B) . (6) 



Here S(b + B) is the S-matrix for a quark-nucleus collision with impact parameter b + B 

— * — * 

where B and b are the impact parameters between the center of gravity of the projectile 
hadron and the center of the nucleus or the quark, respectively. We take the trace over the 
color indices of the quark. The initial density matrix reads, 

6?„(6i, A) = E \°n? J d% d% ... d 2 b n *t (61, b 2 , b n ) b 2 , b n ) . (7) 

n, {in} 

Here we sum over different Fock components of the hadron containing different numbers n 
of (anti) quarks with weight factors \C n \ 2 . We also sum over the initial state polarizations 
and colors of all quarks except for the first one. This makes the matrix Q^ n (bi, b[) diagonal 
in color space. We do not show explicitly the dependence of the hadron light-cone wave 
function on the longitudinal momenta of the quarks, assuming integration over all longi- 
tudinal momenta except for that of the first quark. In the high energy approximation we 
neglect the energy loss of the quark propagating over a finite path in the nuclear medium. 
This effect, if it becomes important, should be treated separately. 

Let us consider the S'-matrix of a quark-nucleus collision in the approximation where all 
the coordinates r*j of the bound nucleons, as well as the intrinsic quark coordinates pj in the 
nucleons, are "frozen" during the interaction time. In this case the S matrix acquires the 
eikonal form, 

S(b + B,fi;{p,/j}i) = E®^ 2 - zi) -&(za ~ z A -\) 

p 

x s 1 {b + B-f Tl ;{p,p} 1 )...s A (b + B-f TA ;{p,p} A ) ■ (8) 

Here we sum over permutations of the nucleons. {p, p}i denotes the set of intrinsic quark 
coordinates and color indices in the i-th nucleon. The single quark-nucleon S'-matrix reads, 



^ (b + B - r Ti \ {p, p]i) = exp 



jj:X a X a U)x(b + B-f Tj ) 
4 3=1 



(9) 



where A a are the Gell-Mann matrices, and index j refers to one of the quarks in the target 
nucleon; 

00 j 2 

X(b) = J%a s (q 2 ).Ub. q ) . (10) 
A 2 q 
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Here Jo is a Bessel function, q is the transverse momentum of the gluon exchanged in the 
t-channel, and A 2 is an infra-red cut off. 

As soon as the initial density matrix Clj n (b, b') is diagonal in the color indices of the quark 
(see above), we can average the product S*S over the colors and coordinates of the quarks 
in the target nucleons in inverse sequence starting from the last one. Then we average over 
the positions of the nucleons {f}, 



+ B, n- {p, fi}i) S(b + b, n- {p, p}, 



E 6(^2 -z x )... <3>(z A - za- 



{P>}/ {r} \ P 

x < s{(V + B- r Tl ) (4(6' + B - f T2 ) . . . (s\(V + B- f T J • s A (b + B - f T _ A 



{p,v}j 



x s 2 (b + B -f T2 )) sx(b + B-f Tl ) 

/{p>}i/ {r - } 



This expression is illustrated in Fig. [I], where we use the two-gluon model for the Pomeron 
for the sake of simplicity. One can see that the Pomerons are enclosed within each other. 




Figure 1: The probability of multiple interactions via one gluon exchange 
for the quark in the nucleus. The dashed line shows the unitarity cut. 



After averaging over the target nucleon wave function including the target quark color 
indices, the functions 

U k (b> + B-f T J+B- f Tk ) = (4(6' + B- f Tk ) ■ s k (b + B - r Tk )) {Mk (12) 
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commute, and the sum over permutations of the nucleon coordinates of the product of 0- 
functions equals one. This expression becomes independent of index k after integration over 
r*fc. The integral can be represented as, 

d 3 r p A {f T , z) U{9 + B - f T ,b + B - r T ) 
- f d 3 rp A (f T ,z)-\ f d 3 rp A (r T ,z) [l - U(b' + B - f T ,b + B - f T ) 



A 



A 



1 " ^T A + B ] a qq (b - b') . 



(13) 



We use the standard approximation of uncorrelated nuclear wave functions, and p A (r) is 
the one-body nuclear density normalized to A. T A (b) is the nuclear thickness function, 



T A {b) = J dzp A {b,z) 



(14) 



and 



*M = 2 I d2r T l-U(b' + B-r T ,b + B-r T ) 



(15) 



is the total cross section for the interaction of a qq dipole of transverse separation p = b — b' 
fllTf with a nucleon. 

Eventually we arrive at the expression for the transverse momentum distribution of a 
quark that has propagated through the nucleus, 



dN q 
d^ 



d bd b' exp ikxib — b') fi* n (6, b') exp 



\^ q (b-b')T A ( b -±^ + B 



■ (16) 



This equation describes the full /cr-distribution of the final quarks, as well as the mean 
transverse momentum squared of an ejectile quark in a hadron-nucleus collision, 

1 



.4 



x exp 



<TB / dzp A (B,z) / d 2 k T K / d%d%n q in (b u b 2 ) 



ik(bi — b 2 ) 



exp 



-0-3,(61 - 6 2 ) T A I bl + b2 + b^ z 



;i7) 



S 



The integration in (|T7|) is easy to perform by replacing k\ exp [i k (bi — b 2 )] by derivatives 
V(6i) ■ V(b 2 ) exp [i k {b\ — b 2 )}. Integrating by parts we arrive at the expression 



1 

A 



(k 2 T )hA = i I d l B I dzp A (B,z) 



+ V(6i) • V(6 2 )exp 



-oo 

1 



Oqq{bi -b 2 )T A (B,z) 



18) 



where we assume that the nuclear radius is much larger than that of the the projectile 
hadron. 

Using a Gaussian for the hadronic wave function, the density matrix of the initial quark 
in the case of a proton beam reads, 



exp 



b\ + b\ 



3 (r 2 ch ) 



(19) 



3vr(r c \) 

where (r 2 h ) is the mean square of the proton charge radius. Note that we do not introduce 
any unreasonably large primordial transverse momentum for the projectile quark, which is 
usually assumed for hard reactions in leading order in the parton model. In the light cone 
approach, the higher order perturbative corrections are already included in the phenomeno- 
logical cross section (|j), and they generate the energy dependence of C(r T ,s). 

In the second term on the r.h.s. of ([18]) one should apply both derivatives to aq q {bi — b 2 ) 
in the exponential. This is because aq q {bi — b 2 ) oc (pi — b 2 ) 2 as b\ — > b 2 , otherwise the 
result is zero. Using Eq. @ for the dipole cross section, one gets from (|T7|) the mean-square 
transverse momentum of a valence quark, 

2 



P A 



3(r c \) 



+ 2C(0,s) (T A ) 



where 



(Ta) 



1 

A 



d 2 bT 2 (b) 



(20) 



(21) 



is the mean nuclear thickness. Thus, with a better derivation we confirm the result of [fL4 |, 
Eq. ®. 
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3 Moliere theory: a more intuitive derivation 



The result we obtained in Eqs. ([jjf), fl20|) might look puzzling. Indeed, in the expansion of 
the exponential in ([Tj|) only the lowest order term oc a qq Ta oc A 1 / 3 survives differentiation 
and contributes to (k^). This observation seems to lead to the conclusion that only single 
rescattering contributes to (k^), while higher multi-fold scatterings do not affect the broad- 
ening of the transverse momentum (compare with p0|). Such an interpretation contradicts 
the intuitive expectation that broadening is caused by multiple interactions in the medium 
resulting in random walk of the particle in the transverse momentum plane. Here we present 



a different, more intuitive derivation for the same result fl20|) based on Moliere's theory [^T 

of multiple interactions in a medium (see also This result shows that broadening is 

indeed caused by multiple interactions rather than by the single rescattering term, which 

has an A-dependence very different from A 1 / 3 . 

We concentrate on nuclear broadening of {kj) and for the sake of clarity neglect the 

primordial transverse momentum distribution of the projectile quark which is responsible 

for the first term on the r.h.s. of (|20|). Evolution of the transverse momentum distribution 

D{k,Ti z ) as function of longitudinal coordinate z is described by the kinetic equation, 

d D(k T , z) ^_ ^ J d2 d a(k T - k' T ) 

a z J d z k T 

Here da jd 2 kx is the differential cross section of quark scattering on a nucleon summed over 

final states and quark colors. The first term on the r.h.s. of (|22]) and the total cross section, 

f ,9, da , . 

otot = J d 2 k T — , (23) 

describe the leakage of quarks from the element of phase space d 2 k^. At the same time, this 
region gains new quarks from other parts of phase space via scattering as is described by 
the second term in (|22|). 

Equation ([22|) is easy to solve switching to coordinate space, 

D(k T ,z) = ^- [ d 2 r T e ir ^^D(r T ,z) . (24) 

2 71 J 
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The solution is (up to the initial condition), 

D(k T ,z) oc ±-J d 2 r T e l ^ kT exp \-(a tot - 7 (r T )) T A (z)] , (25) 



where 



7 (r T ) = Jtfhre^*^, (26) 

7(0) = at* , (27) 

z 

T A (z)= f dz'p A (z r ) . (28) 



We are now in a position to discuss the meaning of expansion of the exponential in (|25|) . 
First of all, the expansion, 

(-l) n 

exp(-a tot T A ) = Kt T A ) n , (29) 

n=0 ^" 

cannot be treated as a multiple scattering series. This is just an attenuation exponential, 
and the particle (quark) cannot be absorbed (knocked out of a phase space cell) twice. This 
simple exponential should not be confused with the Glauber eikonal multiple scattering 
series, which has a nontrivial interpretation in terms of unitarity relation and AGK cutting 
rules IE2. 



At the same time, the expansion of the second term, 7(r^) T A , in the exponent of ( |2~5| ) 
does have the interpretation of a multiple scattering series. Indeed, 



J d 2 rT exp(i %t ■ Ft) exp ^(ry) T4 
S<&) + Ef/f[^^5 (kr - ± k) . (30) 



The n-th term on the sum in the r.h.s. of this equation clearly corresponds to the n-fold 
scattering of the quark. Thus, all multi-fold rescatterings contribute to the shape of the 
/cy-distribution of a quark propagating through a medium. 

Amazingly, the exponent in ( p5|) can be represented as a color dipole cross section, 
making it similar to the results of previous section, 

da 1 

'T 
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*to t - l{r T ) = i / d 2 k T (l - e a ^) = \ a, q {r T ) . (31) 



We have made use of the fact that da/d 2 kx oc l/fcj.. This expression is also derived in more 
detail at the end of this section. 

The appearance of the dipole cross section, as in ([FTP, is the result of an artificial 
construction. Clearly, the object participating in the scattering is not the colored qq dipole 
but rather a single colored quark. Now it turns out that the different parts of this dipole 
cross section have quite different origins. Namely, the first term in the brackets in (|3lD (the 
"1") corresponds to simple attenuation of the projectile quark detected in a given phase 
space cell. However, the second term (exp(ikj- ■ kr)) originates from multiple scattering of 
the quark. 

Thus, if one needs to establish a relation between the expansion of the exponential 
exp[— \(Jqq{rT) Ta] in fllTD and the multiple quark interaction, it would be incorrect to think 
that the n-th order term of this expansion corresponds to the probability W n to have n-fold 
quark multiple scattering, 

W n (r T ) + t±L\ a (r T )T A \ n , (32) 

but, instead, 

W n (r T ) = e~ a ^ TA J^h^Il . (33) 

In contrast to (H), all the terms in (|33|) are positive as they should be. This standard result 
of the multiple scattering theory has been also revived recently in [J23|J . 

Now we are in position to figure out where the mean k\ comes from. We find 

(4> = — / d 2 k T k 2 T [ d 2 r T e rh ^ £ W n (r T ) = k% (n) , (34) 



CTtot 

where 



n=0 



kl = — [d 2 k T e T ^, (35) 
a tot J dkj, 

is the mean transverse momentum squared gained by the quark in a single scattering, and 

(n) = crtot T A (36) 
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is the mean number of interactions of the quark propagating though nuclear thickness T4. 

Equation (|34j) explicitly demonstrates that nuclear broadening of the quark mean trans- 
verse momentum is the result of multiple interactions leading to a random walk in the 
transverse momentum plane. Therefore, it would be incorrect to interpret the result of the 
previous section, Eq. (|T7]), as the contribution of a single scattering, which must contain an 
extra factor exp(— a to t Ta) and has quite a different A-dependence compared to Ta oc A 1 / 3 
of Eq. fll] 



Concluding this section we present a derivation of Eqs. (17) - (|20|) within the model 
of potential scattering. The amplitude for elastic scattering of a particle from a potential 
reads, 



f(k T ) = — I d 2 be^u(b) 

I 7T 



where 



u(b) = 1 - e 



ix(b) 



X(b) 



dzV{b,z) . 



Corresp ondingly, 



a tot = 47rIm/(0) = 2ReJ d 2 buj(b) 
da 



d?kn 



d 2 kn 



da 



J d 2 b (2-2 cosx 



Otot 



(37) 

(38) 
(39) 

(40) 
(41) 
(42) 



as one should have expected for potential scattering. 
Substituting (|l]) into (||) we get, 

j(r T ) = 



d 2 bu*(b + arr) u(b-(l- a)f T ) = a tot - -a qq (r T ) 



(43) 



where < a < 1 is an arbitrary number, and 



2Rey d 2 b jl — exp ix(b + ar T ^j — i\{b — (1 — a)r T ^j j 



(44) 
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is the total cross section for a quark-antiquark pair with transverse separation ry, called 
dipole cross section. In this case a can be interpreted as a share of the total light-cone 
momentum carried by a quark or antiquark. Thus, we confirm Eq. (|31|). 



4 Calculation of the C(0, s) 

As mentioned above, all total cross sections rise with energy, so the factor C(tt,s) does 



also. The energy dependence is steeper towards small ry P^fl and may lead to dramatic 



changes in C(tt, s) compared to the oversimplified estimate C ~ 3 |11. 15]. 

Intuitively, it is clear that a fast quark scatters off the gluon clouds of bound nucleons. 
The gluon density at small Bjorken x relevant to such a high-energy interaction increases 
with 1/x, i.e. with quark energy (see Eq. [|). Said differently, multiple interactions of the 
quark are accompanied by gluon bremsstrahlung. As the energy is increased, more phase 
space becomes available for the radiated gluons. Since gluon radiation enhances transverse 
motion of the quark, (/c T ) should grow with energy (via C(0, s)). 

Thus, the energy dependence of C(tt, s) can be expressed in terms of the unintegrated 
gluon density |TT|| , 



2 



C(0,s) = lJd>k^^F(x,k 2 ). (45) 

where !F(x, k 2 ) = dG(x, k 2 )/d(lnk 2 ) and G(x,Q 2 ) = xg(x,Q 2 ). Here k is the transverse 
momentum that a quark acquires scattering off the gluon cloud of a nucleon. a s (k 2 ) is the 
QCD running constant, which we calculate in the one-loop approximation. 

The value of x for gluons in the r.h.s. of (|45|) is related to the quark-nucleon cm. energy 
squared s. The minimal value of x corresponding to a collinear quark-gluon collision reads, 

Ak 2 

S 

The density of the gluon cloud around the proton should vanish at large distances be- 
cause of the color neutrality of the proton. Therefore, the unintegrated gluon distribution 
T{x, k 2 ) oc k 2 as k 2 — > 0. This provides infra-red stability of the integral in (flop . 
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There are quite a few models for the unintegrated gluon distribution T{x, k 2 ) at high 
k 2 (e.g. see [^]); however, very little information is available in region of small k 2 where 
perturbative QCD cannot be used. Since the integral in ( |45| ) is dominated by soft gluons, one 
should develop a phenomenology for J-(x, k 2 ) and use other observables to restrict possible 
uncertainties. One of the sensitive probes for !F(x, k 2 ) (even more sensitive than C(0, s)) is 
the total hadronic cross section. The simplest case is the pion-proton cross section, which 
is given by the same approach as Eq. (fi5|), 



a^(s) = ^ld 2 k 



a s {k 2 
k A 



1-W H^k 2 ) . 



(47) 



Here F^Jk) 



71 



cxp 



ik(f 2 — fi) 7r) ~ cxp(— k {r 2 )/3) is the two-quark formfactor of 



the pion, and {r 2 h ) n = 0.44 ± 0.01 fm |26[ is the mean-square pion charge radius 



In the Born approximation |27|, [28|, |TT| the gluon density takes a simple form, 

4 



where FP q (k) = (N 



exp 



T(x,k 2 ) 
i k (r* 2 — fi ) 



71 



a s (k 2 ) 1-F£(k) 



(48) 



iVy w exp(— k 2 (r 2 h ) p /2) is the two-quark formfactor 
of the nucleon, and (r 2 h ) p = 0.79 ± 0.03 fm 2 is the mean-square charge radius of the 
proton. In fact, because Coulomb gluons have no partonic interpretation, quark elastic 
scattering through the exchange of two such gluons cannot be expressed in terms of the 
gluon density as given in Eq. (|4]). Therefore, ( |48D cannot be treated as a model for the 
gluon distribution function jF(x, k 2 ) by any means; nevertheless we plot it vs k 2 in Fig. |2] 
(as the dotted curve) to compare with other models. This contribution is independent of x. 
As was mentioned, x-dependence originates from gluon bremsstrahlung. 

Comparing the factor (1 — F£ q (k)) in (f47D to k 2 in (flST), it is clear that the total hadronic 
cross section is much more sensitive to the behavior of JF(x, k 2 ) than C(0, s) is in the region 
of interest, namely at small k 2 . Therefore, reproduction of the total hadronic cross section 
is an important test for any model of the gluon density. At large values k 2 perturbative 
QCD is supposed to be valid, and different models should therefore converge there. 
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k* (GeV*) 

Figure 2: The non-integrated gluon density k 2 ) calculated at the 
quark-nucleon energy s = 500 GeV 2 as function of k 2 . The dashed and 
solid curves correspond to the models GBW and KST. The dotted curve 
shows the two-gluon contribution (j^S). 



Note that one should interpret fl47|) as the inelastic rather than the total cross section 
because it is linear in gluon density, i.e. corresponds to a color octet-octet unitarity cut, 



which does not contain elastic or diffractive states. Therefore, one should add to (f47j ) the 
elastic and single diffraction cross sections to get the total cross section. However, such 
an amplitude linear in the gluon density is subject to unit arizat ion. The first unitarity 
correction, which is quadratic in the gluon density, is just the same elastic and diffractive 
cross sections, but with a negative sign. Hence, this cancels the contributions of the elastic 



and diffractive cross sections, and one can treat (|47|) as the total cross section neglecting the 
higher order unitarity corrections 0(G 3 ), which are known to be quite small at the present 
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energies. 

We will try a few models for the gluon density to calculate the C(0, s), hoping that the 
spread in the results may serve as a measure for the theoretical uncertainty. 



4.1 GBW model 



Golec-Biernat and Wiisthoff |30| suggested a model for the dipole cross section (|2|) which 
saturates at large qq separations and reproduces well the DIS cross sections. 



a qq (r T ,x) = (x 



1 — exp 



T 



(49) 



where R Q (x) = 0.4 fm x (x/x ) x/2 and a = 23.03 mb; A = 0.288; x = 3.04 • 10" 4 . This 
cross section obviously satisfies the behavior suggested by data for F 2 (x, Q 2 ): the smaller 
the separation r^, the steeper the growth with 1/x. However, the value of x is not well 
defined for a dipole of a given energy (there is no problem with the definition of x in our 
next model). 

The gluon density corresponding to this dipole cross section reads J50], 



3 (To 



k R (x)exp 



--Rl(x) k< 



(50) 



16 7r 2 a s (A; 2 ) 

This function is depicted in Fig. ^] by the dashed curve. Note that according to Eq. (|50|) 
the right wing of this curve at large k 2 is a steeply rising function of energy. 

This model faces obvious problems with hadronic cross sections. Indeed, the dipole 
cross section (|^) grows steeply with r T and saturates for r T > Rq(x) at <Tg q (r T ,x) < (To- 



Averaging (f49"D weighted with the pion wave function squared, one never can reach a pion- 
proton total cross section larger than o"o = 23.03 mb. 

In order to calculate the energy dependence of C(0,s) we need to know the value of 
x, which is poorly defined in this model. If we use the minimal value fl4"ED permitted by 
kinematics, the result shown by the dashed curve in Fig. |3] should be an upper bound for 
C(0, s) in this model. We see that C(0, s) is rather steep as a function of energy. 
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Figure 3: Factor C(0, s) as function of the quark energy calculated with 
(fio[). The curve assignments is the same as in Fig. ^, and the thin solid 
curve is the KST curve corrected for gluon shadowing in lead. 



4.2 KST model 



The advantage of the GBW parameterization ([49]) is simplicity and convenience for analyt- 



ical calculations. One may keep this form, but modify it as suggested in ]3TJ to address the 
problems mentioned above. An explicit energy dependence is introduced in the parameter 
do in (f4"9p in a way that guarantees the reproduction of the correct hadronic cross sections, 



00 00 = <Ttrt(s) 1 + 



(51) 



where cx^(s) = 23.6 x (s/so) is the Pomeron part of the irp total cross section P2| , and 
R (s) = 0.88 fm x (s/s )" A/2 with A = 0.28 and s = 1000 GeV 2 is the energy-dependent 
radius. With these parameters the proton structure function F 2 (x,Q 2 ) was calculated [^TJ 
using the non-perturbative distribution functions for the qq component of the photon [|3T 



and the dipole cross section ((49|), ([51]). The results agree well with the data up to Q 2 ~ 
lOGeV 2 , which is sufficient for our interval of k 2 . 

The corresponding gluon density ( |50D is shown as function of k 2 by the dashed curve in 
Fig. with the redefined functions cr (s) and Ro(s). The factor C(0, s) calculated with this 
gluon density is shown by thick solid curve as function of energy in Fig. [| C(0, s) rises with 
energy similar to the prediction of the GBW model. 

Summarizing, the two models under consideration provide quite different values and k 2 - 
dependences for the non-integrated gluon density, as one can see from Fig. |[ Nevertheless, 
in the energy range of the E772/E866 experiment at Fermilab the value C(0, s) ~ 3 — 4 is 
pretty certain and is about twice as big as the simple two-gluon approximation predicts. 

5 Nuclear shadowing of gluons 

An important source of broadening in the transverse momentum of a quark is the gluon 
radiation that accompanies the multiple interactions of the quark in the nucleus. Indeed, 
one can see from Fig. |3| that the factor C(0, s) is enhanced compared to the contribution of 
two-gluon exchange, which corresponds to a quark scattering on Coulomb gluons without 
gluon radiation. 

It is known that radiation in multiple scattering is subject to Landau-Pomeranchuk 



suppression | 33| , which is a coherence phenomenon in radiation. Namely, if the gluons 
radiated due to scattering of a quark off different nucleons are in phase one should add 
up the amplitudes rather than the probabilities. Interferences substantially suppress the 
radiation compared to the classical expectation (Bethe-Heitler approximation). 

The same phenomenon may be treated quite differently as gluon fusion in the parton 
model in the infinite momentum frame of the nucleus, where it is known as nuclear shadowing 
for the gluon density at small x |3Hfl-||38|| . Indeed, only the fast part of parton clouds of bound 



nucleons are squeezed by Lorentz contraction, but the low-x partons may be spread in the 
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longitudinal direction more than the longitudinal size of the nucleus in its infinite momentum 
frame. Therefore, partons originating from different nucleons at the same impact parameter 
overlap and may fuse. This effect diminishes the number of partons at small x. Thus, 
we replace the gluon density in (f45[) by a shadowed one, suppressed by a factor Sa(x, k 2 ) 
compared to the unintegrated gluon density in a free proton, 

C(0, s) = | J d 2 k F(x, k 2 ) S A (x, k 2 ) , (52) 

where Sa(x, k 2 ) is calculated at x — 4k 2 /s. 

The nuclear suppression factor Sa(x, k 2 ) is calculated in |3T| using the light-cone Green 
function approach. At small k 2 , shadowing is controlled by the strong nonperturbative 
interactions of gluons, ensuring that shadowing is nearly independent of k 2 for k 2 < 4GeV 2 . 
This covers the region of k 2 we are interested in, therefore we can safely disregard the k 2 
dependence of Sa{x, k 2 ), except for that which comes from the ^-dependence. The factor 
Sa(x, k 2 ) is depicted in Fig. [| at small k 2 for a few nuclei as calculated in |JT]. We see 
that the onset of gluon shadowing takes place at x < 10~ 2 , which are smaller values of x 
than one needs for the shadowing of quarks (F 2 A (x, Q 2 )). Nevertheless, one can see from 
Fig. |2| that the values of k 2 where the unintegrated gluon density is large are quite small, 
corresponding to very small x where one may expect strong shadowing effects. 

The results of calculation of C(0,s) for lead with the KST model for the unintegrated 
gluon density are shown in Fig. || by the thin solid curve. It turns out that the effect of 
gluon shadowing is not dramatic, i.e. does not exceed 10%. 



6 The effect of final aperture 

The factor C(0, s) for the nuclear broadening of the transverse quark momentum in (|20|) 
appears as a result of integration over kx up to infinity. In reality, the angular acceptance 
of the experimental apparatus restricts the accessible values of kr to be less than some k m . 
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Figure 4: The nuclear shadowing factor Sa(x, k 2 ) for soft gluons for lead, 
copper and carbon (from bottom to top) calculated with the light-cone 
Green function approach J3l|/ . 



Introducing this upper cut off for kr in fll7|) we arrive at a modified value C(s) that has to 
replace C(0, s), 



where 



F(B,r) 



and 



C(s) 



2 7T k m 



oo oo 



dBB drJ x (k m r)F(B,r) 



o o 
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Vq q (r, s) 



9(Tq q (r } S) 

Or 2 



d(Tq q (r, S) 

dr 2 



1 2 



3 { r 1h)v, 
f2(B,r) 



d 2 CTqq(r,s) r2 

Qf2i /y 2 



fi(B,r) 



(53) 



(54) 



h{B,r) = l-(l + a)e~ a 



(55) 
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h(B,r) 



l + a + j e 



(Tg q (r, s) T A {b) . 



(56) 
(57) 



Integration over B and r in (^) needs to be done numerically. The suppression factor 
K{k m ) = C(s)/C(0,s) is plotted in Fig. | as function of k m . This substantially deviates 




0.5 1 1.5 2 2.5 3 3.5 

k m (GeV/c) 

Figure 5: The suppression factor K{k m ) = C(s)/C(0,s) caused by a 
finite aperture calculated with Eq. [5^ as function of k m . 

from unity for a cut off k m < 2 GeV. 



7 Discussion of the results 

We performed calculations for broadening of the transverse momentum of a quark propa- 
gating through a medium based on the light-cone color dipole representation, which treats 
the broadening of transverse momenta as a color filtering of transverse sizes of qq dipoles 
propagating through nuclear matter. It is natural that the broadening is proportional to 
the length of the path times the density of the medium, since multiple interactions cause 
the quark to undergo a random walk in the transverse momentum plane. The main task is 
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to calculate the coefficient C which turns out to be the same as in the dipole cross section 
Cq g ( r T) = C r\ (see Eq. (p0|)). This is not surprising since, in our approach, the k T distribu- 
tion of the final quarks is proven to have the form of an eikonalized dipole cross section. In 
addition to the formal derivation, we present another one based on the Moliere scattering 
theory. This simple and intuitive approach helps to understand how the propagation of a 
single quark can be described in terms of the color dipole cross section. In particular, it 
turns out that a part of the dipole cross section originates from the simple attenuation of the 
quark traveling through the medium, while the rest of the cross section is due to multiple 
interactions of the quark. It is demonstrated that the wide-spread belief that the lowest 
order rescattering leads to a broadening 5{k^) oc A 1 / 3 is incorrect. 

We have also performed a parameter-free evaluation of the effect of broadening. We 
calculated the factor C(0, s) using a phenomenological description of the dipole cross section 
that is adjusted to data for the proton structure function and total hadronic cross sections. 
We found that the effect of /^-broadening steeply rises with energy. 

Several corrections diminishing the effect of broadening, such as nuclear shadowing of 
gluons and the finiteness of the experimental aperture, are evaluated as well. 

We make a few remarks in conclusion: 

• For gluons propagating through a medium, our results for /^-broadening must be 
enlarged by the Casimir factor 9/4. An example of a process relevant to gluon prop- 
agation is heavy quarkonium production. 

• Numerically, our results are in a reasonable agreement with /^-broadening observed 
in production of J/ty and T off nuclei, but are somewhat higher than what is observed 
in Drell-Yan processes 0. A detailed comparison with available data for Drell-Yan 
processes, heavy quarkonium production, dijet and dihadron production off nuclei, 
etc., will be presented elsewhere. 

• Although we concentrate in this paper on the calculation of the nuclear broadening 
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of the mean-square transverse momentum S(k^), we are in a position to calculate the 
full momentum transfer distribution of partons after they escape from the nuclear 
medium. This distribution is given by Eq. (|T6|). 

• Nuclear broadening of the transverse momentum of partons originating from relativis- 
tic heavy ion collisions is enhanced, since the parton experiences multiple interactions 
propagating through both nuclei. On top of that, if the parton is produced at mid 
rapidity with a high transverse momentum it should substantially increase broadening 
relative the original direction if it propagates through a dense matter (i.e. quark-gluon 
plasma) since it probes the density of the medium. This needs measurements with 
two back-to-back hadrons or jets. Thus, the effect of broadening of (k^) can serve as 
a sensitive probe for creation of dense matter in heavy ion collisions. 
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